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Abstract 

We present results of lattice simulations of the Plane Wave Matrix Model 
(PWMM) . The PWMM is a theory of supersymmetric quantum mechanics that 
has a well-defined canonical ensemble. We simulate this theory by applying ra- 
tional hybrid Monte Carlo techniques to a naive lattice action. We examine the 
strong coupling behaviour of the model focussing on the deconfinement transi- 
tion. 



1 Introduction 



The AdS/CFT correspondence [H Ej [3] has emerged as a useful setting for studying 
certain strongly coupled gauge theories. Though this correspondence has been widely 
tested, it is useful to look for examples in which the correspondence can be tested 
directly by studying theories in which tractable calculations can be done directly at 
strong coupling. Lower dimensional examples would seem to provide a useful setting 
for such a study. Recent effort has been devoted to studying the BFSS matrix model 
at strong coupling using both lattice and non- lattice techniques [HOE]. However, 
even that one dimensional quantum theory is still not well defined in the canonical 
ensemble because it has a moduli space The moduli space in this theory is given 
by the eigenvalues of nine commuting matrices that transform under the 5*0(9) R- 
sjTumetry. The existence of these fiat directions means that, at finite temperature, 
the partition function is formally divergent, and it was shown in [7] that when Monte 
Carlo simulations of this theory are performed, this divergence eventually causes the 
simulation to break down. 

The Plane Wave Matrix Model (PWMM) [8] arises as a natural next choice for sim- 
ulation because it includes supersymmetry preserving mass terms that lift the moduli 
space, and give a discrete set of vacua. The theory still has enough supersymmetry 
that it has a known gravity dual at zero temperature [5]. Symmetry reduces the dual 
supergravity problem to a two dimensional one, and the geometries dual to the discrete 
set of field theory vacua can be put in one-to-one correspondence with a set of axisym- 
metric electrostatics configurations, the potential of which determines the supergravity 
fields, and was solved in [10]. Moreover, this theory has a number of large limits 
that give theories of physical interest: M-theory on a plane-wave [8], little string theory 
on [n], and TV = 4 SYM on RxS^ [I2]. Studying the PWMM at strong coupling, 
therefore, is a first step toward direct simulation of these theories. 

The PWMM is also a good candidate for establishing benchmarks for studies of 
lattice supersymmetry in general. It is simple enough to allow for direct study using 
a naive lattice action, is well-defined in the canonical ensemble and of course being 
one dimensional is numerically very tractable. As such it is complementary to studies 
of supersymmetric theories in higher dimensions using new lattice formulations which 
retain some exact supersymmetry at non zero lattice spacing - see the recent review 
[T3] . For example it should be possible to cross check the results of direct simulations 
of A/" = 4 super Yang-Mills in certain limits with results derived from simulations of 
naive lattice discretizations of the PWMM. 

In this paper we will present the results of simulations of the PWMM on the lattice. 
We use a naive lattice action, and concentrate on studying the Hagedorn/deconfinement 
transition in the model. We study the theory at fixed temperature, measured in units 
of the mass deformation, as a function of the 't Hooft coupling (measured in the same 
units). We used a quenched approximation to explore the dependence of the critical be- 
haviour on the rank of the gauge group and the number of lattice points. The quenched 
approximation is much less computationally demanding and can be used to estimate 
reasonable choices of parameters for simulating the full theory. We establish that we 
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can get a reasonable approximation of the continuum large behaviour with modest 
values of both TV and the size of the lattice. We then simulate the full theory. The 
Hagedorn/deconfinement transition in the PWMM has also been studied previously at 
weak coupling 't Hooft coupling in the large limit, and we compare our results to 
extrapolations of those from weak coupling. 

The remainder of this paper is organized as follows. In section [2] we review the 
PWMM, its vacuum solutions, and various limits. In section [3] we present our lattice 
action and discuss the simulation method. Section H] contains our main results. We 
review the weak coupling results and contrast them with the results of our simulations. 
In section [5] we discuss possible completions of the phase diagram. 



2 The PWMM 

The Plane- Wave Matrix Model (PWMM) is a gauged matrix quantum mechanics the- 
ory with sixteen supercharges. It has an action that can be written as 



i i<j i 



in which X^ are nine scalar matrices that sit in the adjoint representation of SU{N), 
i.e. X* = Yla=i^ ^i'T"', where are anti-Hermitian generators that are normalized 
to Tr T"T^ = —6ab', ^1/ is a Majorana fermion that is also in the adjoint representation 
of the gauge group; Tr is over the gauge indices; D is a gauge covariant derivative; A 
is the 't Hooft coupling, given by A = g'^N; and 

7123 = ^(^ijkliljlk ■ (2) 



The terms in the first line of ([T]) are those in the BFSS matrix model [Hj. The 
terms in the second line give supersymmetry preserving masses to the fields and add 
a Myers term. The addition of these terms gives this model two technical advantages 
compared with BFSS. Firstly, in the limit yU — )■ oo the model becomes weakly coupled 
and can be studied perturbatively [Hj. Also, the addition of these terms lifts the 
moduli space of the BFSS model giving instead a discrete set of vacua [8]. This is 
important for numerical simulation because it means the theory is well-defined in the 
canonical ensemble. If we consider the scalar fields X* for i = 1, 2, 3 we can write the 
potential as 

(/iX^ + i=6,,,[X^X^l)2 , (3) 

which is minimized when the scalar fields X* are proportional to generators of SU(2). 
The discrete vacua correspond to the various ways of forming N x N matrices that 



V = -Tt 
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generate SU(2), and can be put in one-to-one correspondence with partitions of into 
a sum of natural numbers. It was shown in [16] that the vacua sit in certain doubly 
atypical representations of the SU{2\4) supersymmetry algebra, which protect them 
from quantum corrections, and thus persist at strong coupling. 

The identification between the field theory vacuum states and regular type IIA 
supergravity solutions with S'[/(2|4) supersymmetry was found in [5]. SU{2\4:) contains 
a bosonic R x SO (3) x SO (6) that reduces the type IIA supergravity problem to a two- 
dimensional one. Imposing supersymmetry and regularity reduces the two-dimensional 
problem to solving a two-dimensional Laplace equation. The vacua of the PWMM are 
encoded by choosing appropriate boundary conditions. Using the fact that the potential 
in electrostatics obeys the Laplace equation, it was shown [9] that regular supergravity 
solutions correspond to axisymmetric configurations of charged conducting discs. This 
class of electrostatics problems was solved in [TO] . 

The PWMM has a number of large A^ limits in which it describes interesting physics. 
See [TT] for a summary; we will mention a few here hereQ 

1. If we consider a limit in which 

N ^oo ^ fixed, (4) 

this model should describe M-theory in a maximally supersymmetric plane-wave 
background. This is because the matrix theory conjecture suggests this model 
describes the DLCQ of M-theory on the plane wave with compact null momentum 
given by 

,2 , A^ A^ 
f^lpP = — oc — , (5) 

where R is the size of the null circle. The limit (jl]) decompactifies the compact 
null circle. 

2. The 't Hooft limit, in which A^ — )■ oo with A fixed, allows for the pertur- 
bative study of the theory if A is small. It is in this limit that the Hage- 
dorn/deconfinement transition can be studied analytically, as we will review be- 
low. 



3. It was shown in [TT] that if we consider expanding around the trivial vacuum 
solution and take a limit in which, asymptotically, 

In^ A^ 

N ^oo —— fixed, (6) 
A 

the PWMM describes type IIA little string theory on S^. Little string theory 
is the theory that describes the degrees of freedom on NS5-branes, and is of 
particular interest because it exhibits some of the features of string theories (e.g. 
T-duality), but not others (e.g. gravity). 



-•^For links between the PWMM and other field theories see [T71 118) . 
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4. Using arguments from the gravity side [9] or from the field theory side [12] TV = 4 
SYM on i? X S"^ can be reahzed by expanding about a particular vacuum of the 
PWMM. (For a full discussion, see also [19], [20].) Expanding about a vacuum 
consisting of representations of SU{2) of dimension n+j, for some integers n, and 
j = 1, 2, 3, . . . , m, each with k copies. If the PWMM has gauge group SU{N), 
then we must have 

N = -km{2n + m + l). (7) 
Taking a limit — )■ oo in which 

m 

n — )■ oo m — )■ oo > 0, (8) 

n 

gives A/" = 4 SYM on R x with gauge group SU{k). 

It is in principle possible to extract the physics of these models by studying the 
behaviour of the model along appropriate curves in the X — N plane, and extrapolating 
to large A^. In this paper we seek to establish benchmarks for this procedure by 
performing a lattice simulation of the PWMM at strong coupling and finite A^. 

3 Lattice Action 

To find our lattice action, we start with the PWMM action ([1]). The bare action has 
three parameters with unit mass dimension As, f3~^ and The gauge field and scalars 
also have unit mass dimension, and the fermions have mass dimension |. We would 
like to make the following rescaling: 

X ^ fiTX D^fiTD t^— ^^(/iT)i^, (9) 

where T is a dimensionless number that we will shortly take to be the number of lattice 
points. We can then write the action as 



ai iri - 2_^{uy^i)- - 

i<j 

9 

i=l i=4 

(10) 



Note that all of our fields have been rendered dimensionless. After integrating over the 
fermions we have the following partition function 

dAdXPi{0)e-^^, (11) 
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with 



. (12) 



2^2 ^ 42^2 

1=1 i=4 



and 



= J'D + Y,l^[X\■] + ^ll23. (13) 



We discretize the theory on a lattice with spacing a, and replace the J dt with a J2f 
We use the derivative operator 

where 

= l(f/(t)^(t + a)f/t _ ^(t)). (15) 
a 

D~ is the adjoint of D^. The lattice action is 

T-l 



and 



m=0 i i<j ^^g^ 

i=l i=4 

O = 7°D^„ + Y,li[XL -limn + ^7l23lmn. (17) 

i 

The overall dimensionless lattice coupling is then 

- = ^^^ (18) 

where /3 is the inverse temperature, and T is the number of lattice points. The contin- 
uum limit is defined by taking T — >■ oo with ^^^^ fixed. 

We simulate this model using Rational Hybrid Monte Carlo (RHMC) [2T] . Since the 
fermions appeared quadratically, we were able to formally integrate them out yielding 
the Pfaffian of the operator O. Using the fact that Pf {O) = det O, the Pfaffian can 
be computed by introducing a path integral over pseudofermions F 

PfO = y"[Z}F]e-^'(^'^)"^^. (19) 
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It is still challenging to compute ((9^(9)" 4, but it can be done efficiently by making 
a rational approximation, and making use of iterative Krylov subspace methods to 
solve the resulting set of shifted linear systems [22]. Fictitious conjugate momenta 
are introduced for each of the variables, and the resulting system is simulated by us- 
ing a combination of molecular dynamics integrations of the equations of motion and 
Metropolis Monte Carlo [231 IM]- The main computational expense derives from simu- 
lating the fermionic Pfaffian, and so it will be useful at times to consider the quenched 
approximation (without fermions). This will allow us to estimate what happens at 
much larger values of than we can when the fermions are included with equivalent 
computational effort. 

4 Hagedorn/Deconfinement Transition 

To get an idea of how well our finite N simulations work we will compute a simple 
observable in the matrix model: the Polyakov loop. There are good reasons to compute 
it. 

• It is a simple gauge invariant observable that is straightforward to define on the 
lattice. 

• It is an order parameter for deconfinement, and is expected to signify a geometric 
transition in the dual gravity theory. 

• It has been computed in the PWMM at weak coupling. 

In lattice gauge theory, holonomies of the gauge field are natural observables. In our 
one-dimensional theory, the holonomy to compute is the one around the thermal circle. 
It can serve as an order parameter for deconfinement because it measures the energy 
cost of adding an infinitely massive external charged particle. In the confined phase, 
the trace of the Polyakov loop has vanishing expectation value. This is because it is a 
unitary matrix, and in the confined phase its eigenvalues are uniformly distributed on 
the unit circle. In the deconfined phase this symmetry is broken, the eigenvalues begin 
to clump and the Polyakov loop develops a vev. 

It is also interesting from the point of view of gauge/gravity duality. In gravity, 
there is a well-known phase transition between thermal AdS space and asymptotically 
AdS-Schwarzschild black holes [25]. Black hole thermodynamics implies that this phase 
transition comes with an increase in entropy. It was argued [26] that this fact, along 
with the fact that the Wilson line computed via gauge/gravity duality ceases to obey 
an area law in the high temperature phase, means this bulk gravity transition corre- 
sponds to deconfinement in the gauge theory. The increase in gravitational entropy 
also matches nicely with this picture because in the confined phase the free energy is 
expected to be of order one, whereas in the deconfined phase it is expected to be of 
order A^^. For a discussion of general aspects of this, see [27]. We therefore expect that 
by studying the deconfinement transition at strong coupling in the matrix model we 
are capturing a geometric transition in the bulk dual. 
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Figure 1: Schematic picture of the phase diagram of the PWMM at weak couphng. For 
above /i/?/^ (at low temperature) the theory is confined. As is lowered through 
uPh the theory undergoes a Hagedorn/deconfinement transition. The transition tem- 
perature can be computed at weak coupling; the calculation in the quenched theory is 
reviewed in section I4.1.H and section 14.2.11 presents the result in the full theory. 



In [271 [28] it was shown that in certain large gauge theories, including the 
PWMM, it is possible to have deconfinement transitions at weak coupling. Let us 
discuss briefly why this occurs, using an argument from [27] . Consider a theory of two 
bosonic matrix oscillators with unit energy in the adjoint representation of a gauge 
group with rank A^. Gauge invariant states are formed by the action of traces of 
matrix creation operators on the vacuum. Let us consider single trace states, which 
dominate in the large A^ limit. The number of single trace states with energy E is 
bounded above by 2^, because there are E positions in the trace that can be filled 
by either oscillator. Some of these states are related by the cyclicity of the trace, 
but there are at most E — 1 such states, and so the number of states with energy E 
is bounded below by This lower bound is sufficient to produce an exponentially 
growing density of states, which indicates that there will be a Hagedorn transition. 

The existence of results at weak coupling suggests where we should look for the 
phase transition at strong coupling. A schematic picture of the phase diagram at weak 
coupling is shown in figure [TJ We will see below that the weak coupling analysis suggests 
that fif3H decreases as the coupling is turned on. As a result, for our simulations, we 
will pick some fixed < fifSn and look for the transition by picking various values 
of the coupling constant. We will first examine what happens in the quenched theory, 
and thereafter present results for the full theory. 
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4.1 Quenched Theory 



Here we will first present results for the quenched theory at weak coupling, which can 
be culled from [281 129], and then the results of our simulations. 



4.1.1 Weak Coupling 

By counting gauge invariant states, it was argued in [28] that a matrix model with 
action 



(20) 



J-- 

with bosonic matrices Xj in the adjoint of SU{N), and A being an SU{N) gauge field, 
has a partition function that, in the large limit, is approximately 

This partition function clearly diverges when the denominator vanishes, and it was 
argued in [2S] that the temperature where this happens is the Hagedorn temperature. 



g-fe^j = 1. (22) 



=1 



If we consider the various fields in the PWMM and decompose them into repre- 
sentations of SU{2\A) which contains a bosonic SU{2) x SU{4), there are scalars that 
transform as a singlet under SU{A) and in the 3 of SU{2) and scalars that transform as 
a singlet under SU{2) and in the 6 of SU{4:). If we expand the quenched theory about 
the trivial vacuum state, we then have three matrices with u = fi and six matrices 
with u = fi/2. Writing xh = e~^"^^'^ we find the Hagedorn temperature is given by 
the solution of 

3xjj + 6xH = 1, (23) 
where the admissible solution satisfies < xh < 1- This gives xh ~ + that 

^13 H = -2 log ( -1 + ^ ) ^ 3.732528074 . (24) 



Since the only massless field in ( 120|) is the holonomy of the gauge field around the 
thermal circle, it is possible to integrate out the matrices to find an effective action. 
This effective action is just the effective action for the Polyakov loop. At one-loop 
level it was shown in [28] that there is a deconfinement transition, whose temperature 
coincides with the Hagedorn temperature. 

It is also possible to calculate the effective action of the Polyakov loop to higher 
order in the coupling constant. This has been carried out in the full theory to two-loop 
order in [29], and three- loop order in [30]. The two-loop calculation is sufficient to 
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determine the correction to the critical temperature, whereas three-loops are needed 
to determine the order of the phase transition. 

The next order correction to the Hagedorn temperature comes from the planar two- 
loop effective action for the Polyakov loop. In this case, using the results of [25], the 
partition function gets corrected tcH 

(25) 



1 - - 6y2 - Xg(y) logy 

with 

giy) = A8y\y^ + 3y' + 7y^ + 13), (26) 

where y = e~^^^^, and we write a dimensionless 't Hooft coupling A = A//i^. Let us 
define 

z{y) = 3y^ + (27) 
The correction to the transition temperature is given by 

yH = yPil-c\nyP~X), (28) 

or, equivalently, 

/i/3H = -41og(|/g))(l-cA), (29) 
for a constant c that is determined by 

. 9{yP) 



(30) 



From ( l29l) it can be seen that c > decreases the inverse temperature as the coupling 
is turned on, whereas c < increases it. Evaluating this at = ^/xH gives 



c 49.04614622 , (31) 

which means that the inverse temperature at the transition decreases as the coupling 
is turned on. 

4.1.2 Strong Coupling 

Here we will study the quenched theory by expressing all length scales in terms of 
/X. The quenched theory is useful for exploring some of the key issues of our lattice 
methods. The PWMM describes interesting physics in various limits in which the rank 
of the gauge group, A^, is taken to infinity. Since the phase space of the model grows 
very quickly with it is necessary to understand if interesting physics is within the 
range of that is computationally accessible. Moreover, we are also limited by the 
lattice approximation itself and we would like to determine how many lattice points 
are required to obtain a reasonable approximation to the continuum. 



^There is also a correction of order A coming from the non-planar diagrams, but it does not 
contribute to the correction to the Hagedorn temperature, so we don't include it. 
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Figure 2: Tr e*'^'^|) in the quenched theory for gauge group SU{N) with various N, 
and T = 5 lattice points. The results agree quite well for all choices of above the 
transition, and well below the transition. The greatest difference is near the transition, 
where it appears to be sharper in the case of larger A^, as expected. 



Figure [2] shows a plot of the Polyakov loop in the quenched theory as a function of 
the lattice coupling measured in units of n, with the inverse temperature = 1, for 
various with five lattice points. The plot shows clear evidence of a deconfinement 
transition. Moreover, it can be seen from this plot that the critical value of the coupling 
does not seem to be very dependent on N. This is very suggestive that the critical 
coupling in the large theory can be reasonably approximated by simulations with 
modest values of A^. Note, however, the transition does appear to be much sharper for 
larger A^, as expected. 

In figure |3] we plot the Polyakov line in the quenched theory as a function of the 
coupling at fixed = 1 and N = 5 for various lattice sizes, T. The plot indicates 
that the critical value of the coupling does not depend strongly on the lattice size, and 
that T = 5 provides a reasonably good approximation to the continuum. 

Let us compare these results with those at weak coupling. If we extrapolate the 
weak coupling results to = 1, we can use fl29|) to estimate the critical value of the 
coupling. We find log(A//i^) ~ —1.826043339. Comparing this value to the plot in 
figure [21 we see that the critical value of the coupling differs considerably. 

There are various possible reasons for this difference. One is that the weak coupling 
computation was done in the 't Hooft limit expanding about the trivial vacuum state. 
The results of our simulations, in contrast, should include the effects of all of the various 
vacua, and are for finite A^. Thinking about our results from the weak coupling point 
of view, as we change A^ our results could be sensitive to both the number of oscillators 
above each vacuum state, and the number of vacuum states, which change with A^. 
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Figure 3: (| Tr e*^"^|) in the quenched theory for various T, N = 5. 
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P|) = 




Figure 4: Schematic picture of the phase diagram from our lattice results. We study 
the theory along the blue dotted line, and identify a critical coupling, Ac. At strong 
coupling (at this temperature) the theory is confined. It undergoes a deconfinement 
transition at Ac. Compare this with the phase diagram at weak coupling in figure [TJ 



11 



We would expect that both of these effects would show up in the dependence of our 
results, yet figure [2] shows the only prominent N dependence is in the sharpness of the 
transition. 



4.2 Full Theory 
4.2.1 Weak Coupling 

The thermodynamics of multi-matrix models with fermions was also considered in 
It was shown that for free gauged matrix quantum mechanics with a collection of 
bosonic and fermionic matrices with frequencies Uj the partition function in the large 
limit is approximately 

^ -—4 ^- (32) 



Ea 



e 



-/3u} 



I.e. fl2T]) continues to hold, however the sum now runs over the frequencies of both the 
bosonic and fermionic matrices. Again, the inverse Hagedorn temperature is given by 
the place where the partition function diverges. Moreover, as in the quenched case this 
can be shown to coincide with the deconfinement transition temperature computed 
from the one loop effective action for the Polyakov loop. 

If we again consider classifying the fields by their transformations under SU (2|4), in 
addition to the bosonic scalars we had in the quenched case we now also have fermions 
that transform in the fundamental representation of both SU(2) and SU{A). If we 
expand about the trivial vacuum state, then, we have, in addition to the three bosonic 
matrices with uj = n and six with u = fi/2, eight fermionic matrices with u = 3/i/4. 
Writing y = e~^^^^ we find the Hagedorn temperature is given by the solution of 

^Vh + 8yl + Qyjj-l = {Svh - i){yH + if = 0, (33) 

where y = yn the transition and the admissible solution satisfies 1. This 

gives i/ff = |, so that 

fx(3H = 4 log 3 ^ 4.394449156 . (34) 



Comparison with the result in the quenched theory (1241) gives a difference of about 
15%. 

As in the quenched case, the next order correction to the Hagedorn temperature 
comes from the planar two-loop effective action for the Polyakov loop. Here, again 
using results of [29], the partition function gets corrected to 

(35) 



1 - ziy) - Xgiy) logy 
with 

z{y) = 3y^ + 8y' + Gy' g{y) = A8y\l + y)\l + y^) . (36) 

As in the quenched case, the correction to the Hagedorn temperature, written in the 
form ([29]) gives ([30]) with 

c = ^ • (37) 
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Figure 5: {|| Tr e*'^'^|) In the full theory with dynamical fermions. The gauge group is 
SU{N) for = 3, 5, and we are using five lattice points. 



Comparing this to the result in the quenched theory, ( l3T]) . shows that the correction to 
the inverse temperature of the transition is stronger by about an order of magnitude 
in the quenched case compared to the full theory. 

4.2.2 Strong Coupling 

We study the PWMM at strong coupling using lattice techniques, as in the quenched 
case. However, the inclusion of the fermionic Pfaffian substantially increases the com- 
putational expense. Fortunately, the results of the quenched simulations suggest that 
five lattice points provide a reasonable approximation to the continuum, at the tem- 
perature we are considering. 

Let us compare these results with those at weak coupling. If we extrapolate the weak 
coupling results to = 1, we can use ( l29l) to estimate the critical value of the coupling. 
We find log(A//i'^) ~ —0.5154039782. Comparing this value to the plot in figure |5l 
we see that the actual critical value of the coupling differs considerably. Another 
interesting order parameter for the transition is shown in figure |6] which plots the 
expectation value of the energy in dimensionless units as a function of the dimensionless 
't Hooft coupling. This is essentially given by the value of the bosonic action since 
the contributions of the fermions can be computed exactly by a scaling argument. 
Clearly, for strong coupling or small mass the vacuum is supersymmetric = 
while supersymmetry apparently spontaneously breaks for small 't Hooft coupling. 

If we compare the critical coupling of the full theory as plotted in figure |5] with 
that of the quenched theory plotted in figure [2] we see that there is little difference. 
This might be surprising, considering that extrapolation of the weak coupling results 
suggested they should differ by more than an order of magnitude. We will comment 
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Figure 6: Energy ^ vs logX/fi^ for SU{5) 



on this in the next section. 

A potential pitfall of our numerical simulations is the notorious sign problem. In 
figure [7] we plot the Pfaffian phase for the theory with gauge group SU{3). The 
expectation value of the phase is very close to zero, as expected, but more importantly, 
the uncertainty in the phase is very small. 



5 Discussion 

We have presented results on lattice simulations of the PWMM. We argued that the 
mass terms that deform the BFSS matrix model to the PWMM render it particularly 
suitable for lattice study by lifting the moduli space and yielding a discrete set of 
vacua. We have focussed on studying the Hagedorn/deconfinement transition in the 
model. By simulating the model at fixed temperature over a range of coupling, we have 
shown that the model exhibits a deconfinement transition when the 't Hooft coupling 
is of order one. We found that critical value of the coupling in the quenched case 
was of the same order as in the full theory. At first glance this might be surprising 
because naively extrapolating the weak coupling result to the temperature at which 
we did our simulations would suggest a difference of more than an order of magnitude. 
However, at sufficiently high temperatures we expect the dynamics to be dominated 
by bosonic zero modes on the thermal circle [31], in which case the critical coupling 
of the full theory should converge to that of the quenched theory. The fact that our 
results seem to be indicate this convergence is setting in at = 1 suggests that in 
future investigation of the phase diagram, it is sufficient to study the quenched theory 
if we are interested in inverse temperatures that sit below (3fi = 1 in figure [HI 

In principle, the thermal ensemble should sample states near each of the discrete 
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Figure 7: The phase of the Pfaffian for the full theory with gauge group SU{3), fixed 
= 1, and five lattice points. The phase is clearly small in the region of parameter 
space we are exploring. 

vacua. In the large limit, since the number of vacua grows exponentially with A^, we 
expect that the moduli space problem encountered in simulations of the BFSS model 
will reappear in the PWMM. The effect of the additional vacua can be enhanced or 
suppressed depending on the physical question of interest using umbrella sampling. 
Though we focussed on the deconfinement transition, and the Polyakov loop is degen- 
erate for non-trivial vacua, we did not seem to need umbrella sampling to restrict to 
the trivial vacuum. A full study of the behaviour of the scalar fields is interesting, 
however, and we leave this for future work. 

Though our simulation does not cover the entire /3 — A plane, it is tempting to 
speculate that the phase diagram is of the form in figure El This seems most natural 
from the point of view of gauge/gravity duality. If the transition at strong coupling was 
not connected to that at weak coupling, for example, that would suggest the existence 
of some new gravity solution. 

We have focussed on one particular aspect of the finite temperature strong coupling 
dynamics of the PWMM, however, there are many other questions of interest; we will 
mention two. In the limit of yU — at fixed f3 and A, the PWMM becomes the 
BFSS matrix model. This model has a well-defined dual black hole solution, and it 
would be interesting to explore this limit to determine if the anticipated black hole 
thermodynamics is reproduced by the strongly coupled matrix model. Also of interest 
is to simulate the theory in the limit in which the PWMM becomes AT = 4 SYM on 
R X S^; results of non-lattice techniques in this direction were presented in [32]. 
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(1^1) = 
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Figure 8: Hypothesized schematic picture of the phase diagram of the PWMM. Our 
simulation scans along the blue dotted line where we are able to identify a deconfine- 
ment transition. The heavy solid black line represents the known weak coupling result. 
The dashed black line represents a suggestion of what the phase diagram could look 
like interpolating between these points. 
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